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New two dimensional systems like surface of topological insulator and graphene offer a possibil¬ 
ity to experimentally investigate situations considered "exotic” just a decade ago. One of those is 
the quantum phase transition of the "chiral” type in electronic systems with relativistic spectrum. 
Phonon mediated ("conventional”) pairing in the Dirac semimetal appearing on the surface of topo¬ 
logical insulator leads to transition into a chiral superconducting state, while exciton condensation 
in these gapless systems has been envisioned long time ago in the physics of the narrow band semi¬ 
conductors. Starting from the microscopic Dirac Hamiltonian with local attraction or repulsion, the 
BCS type gaussian approximation is developed in the framework of functional integrals. It is shown 
that due to an "ultra-relativistic” dispersion relation there is a quantum critical point governing 
the zero temperature transition to a superconducting or the exciton condensed state. The quantum 
transitions that have critical exponents very different from the conventional ones. They belong to 
the chiral universality class. We discuss the application of these results to recent experiments in 
which surface superconductivity was found in topological insulators and estimate feasibility of the 
phonon pairing. 

PACS numbers: PACS: 74.20.Fg, 74.90.+n, 74.20.Op 


I. INTRODUCTION 

Topological insulator (TI) is a novel state of matter 
in materials with strong spin - orbit interactions that 
create topologically protected surface states 1 . The elec¬ 
trons (holes) in these states have a linear dispersion re¬ 
lation, see Fig. 1, and can be described approximately 
by a (pseudo) relativistic two dimensional (2D) Hamil¬ 
tonian. The system realizes an ’’ultra-relativistic” 2D 
electron or hole conducting liquid along with much bet¬ 
ter studied graphene^, a 2D one layers sheet of carbon 
atoms that became a paradigms example of the Dirac 
semi-metal. In the context of graphene certain quan¬ 
tum phase transitions were theoretically contemplated. 
The superconductivity in graphene has been repeatedly 
considered-; however despite great experimental efforts 
was never achieved. The 2D dimensional phonons seem 
to be unable to overcome strong Coulomb repulsion in or¬ 
der to create a Cooper pair. The same can be said about 
attempts to achi eve exciton condensate in graphene that 
was proposed 4 ® even before its discovery. Apparently 
the repulsion is not strong enough either to create stable 
electron - hole bound stated. The surface of topological 
insulator therefore became a prime candidate to realize 
the quantum transitions. 

It is known for a long time that similar 2D and 
quasi-2D metallic systems like the surface metal on twin 
planed^, layered materials (strongly anisotropic high T c 
cuprateiP or organic superconductor^ may develop 2 D 
(surface) superconductivity. This phenomenon became 
known as ’’localized superconductivity”®!. Since best 


studied TIs possess a qui te sta ndard phonon spectruid^, 
it was predicted recently 12 ^ that they become super¬ 
conducting. The predicted critical temperature of order 
of IK is rather low (despite a fortunate suppression of 
the Coulomb repulsion due to a large dielectric constant 
£ ~ 50), the nature of the ’’normal” state (so-called 2D 
Weyl semi-metal) might make the superconducting prop¬ 
erties of the system unusual. The ultra-relativistic nature 
manifests itself mostly when the Weyl cone is very close 
to the Fermi surface. Especially interesting is the case 
(that actually was originally predicted for the [ 111 ] sur¬ 
face of Bi 2 Te 3 and Bi^Se-.p^) when the chemical poten¬ 
tial coincides with the Dirac point. Although subsequent 
ARPES experiments 1 show the location of the cone of 
surface states order tenths of eV off the Fermi surface; 
there are experimental means to shift the chemical po¬ 
tential, for example by the bias voltag^®. 

Unlike the more customary poor 2D metals with sev¬ 
eral small pockets of electrons/holes on the Fermi sur¬ 
face (in semiconductor systems or even some high T c 
material^), the electron gas TI has two peculiarities es¬ 
pecially important when pairing is contemplated. The 
first is the bipolar nature of the Dirac spectrum: there 
is no energy gap between the upper and lower cones. 
The second is that the spin degree of freedom is a ma¬ 
jor player in the quasiparticle dynamics. This degree 
of freedom determines the pairing channel. The pairing 
channel problem was studied theoretically on the level of 
the Bogoliubov-deGennes equation 16 . Both s-wave and 
p -wave are possible and compete due to the breaking of 
the bulk inversion symmetry by the surface. Various pair- 
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ing interactions were considered to calculate the DOS 
measured in Cu x Bi 2 Se 3 using self-consistent analysis®. 
As mentioned above the most intriguing case is that of 
the small chemical potential that has not been addressed 
microscopically. It turns out that it is governed by a 
quantum critical point (QCP) LS . 

The concept of QCP at zero temperature and varying 
doping constitutes a very useful language for describing 
the microscopic origin of superconductivity in high T c 
cuprates and other ’’unconventional” superconductortP. 
Superconducting transitions generally belong to the U (1) 
class of second order phase transition}®, however it was 
pointed out a long time agcPSl that, if the normal state 
dispersion relation is ’’ultra-relativistic”, the transition at 
zero temperature as function of parameters like the pair¬ 
ing interaction strength is qualitatively distinct and be¬ 
longs to chiral universality classes classified in ref.®} The 
term ’’chiral” appears following the corresponding discus¬ 
sion of the well studied both theoretically and experimen¬ 
tally chiral symmetry breaking transition in Quantum 
Chromodynamics. Attempts to experimentally identify 
second order transitions governed by QCP in condensed 
matter included quantum magnets^, superconductor - 
insulator transition }- 22 i and more recently exciton con¬ 
densate in graphen e! 2 ! 23 ! and other Dirac semi - metals 
including TI. In the last two cases the broken symmetry 
is also often termed ” chiral”. 

Exciton condensation is a very old concept in low di¬ 
mensional narrow gap semiconductor physics. The best- 
studied exciton condensate is the quantum Hall bilayer 
at half-filled Landau levels^. Here ingenious methods 
had been developed to separately contact the two layers 
so that one can directly probe the order parameter via 
counterflow superfluidity along the layers and tunneling 
between the layers. The same idea was extended to bi¬ 
layer graphene and recently to Tl 25 . 

In addition Dirac semimetal was realized in cold atom 
system®! (following the realization in 2D known as 
the ’’synthetic graphene”). Interestingly the sign and 
strength of the interaction can be controlled. The Dirac 
semimetal in optically trapped cold atom system}® is 
well suited to study this fascinating phenomenon. The 
Dirac semimetal in optically trapped cold atom}®! of¬ 
fers a well controllable system in which this phenomenon 
occurs both for repulsive interaction (chiral symmetry 
breaking) and the attractive one (superconductivity). 

In this paper the quantum phase transitions in Dirac 
semimetal due to local interactions both attractive (su¬ 
perconductivity) and repulsive (exciton condensation) 
are studied with emphasis on their distinct critical¬ 
ity. The critical exponents belong to chiral universality 
classes that are identified. 

In Section II the general framework that allows to 
study the surface superconductivity and exciton conden¬ 
sation in general Dirac semi-metal (TI, not necessarily 
time reversal and reflection invariant or some other of 
the numerous systems being identified recently) with a 
general local interaction is presented. The local coupling 


strength g , chemical potential p, and temperature T will 
be kept general (g is negative for repulsion leading to the 
exciton condensation or positive leading to superconduc¬ 
tivity). Since the symmetry analysis is crucial, we first 
discuss the space and spin rotations. In Section III we 
concentrate on the simplest time reversal and reflection 
invariant Dirac model and identify its spontaneous sym¬ 
metry breaking patterns. In Section IV the phase di¬ 
agram for g > 0 is obtained for arbitrary temperature 
T and chemical potential much smaller than the Debye 
energy Tp. The latter condition is the main difference 
from the conventional BCS model in which g » Tp. A 
quantum critical point at T = g = 0 when the coupling 
strength g reaches a critical value g c dependent on the 
cutoff parameter Tp. We concentrate on properties of 
the superconducting state in a part of the phase diagram 
that is dominated by the QCP. Various critical exponents 
are obtained. In particular, the coupling strength de¬ 
pendence of the coherence length is £ oc (g — g c ) v with 
v = 1 , the order parameter scales as A a (g — g c )^, 
fj = 1. For the repulsion similar transition occurs in 
the exciton channel in Section V. The critical exponents 
beyond mean held and experimental feasibility of super¬ 
conductivity are discussed in Section VI. 


II. GENERALIZED MEAN FIELD 
APPROXIMATION FOR LOCAL FOUR - FERMI 
INTERACTIONS 

A. Hamiltonian and the partition function 

We consider the second quantized electron Hamilto¬ 
nian via four-Fermi local coupling of strength g 

H = J d 2 rip ] a (r) K a pipp (r) -|v4, (0 tpp (r) tpp (r) ip a (r); 

K a p (V) = £ a p (V) — fj,5 a p, (1) 

where space is two dimensional, r = {x, y} and /i is the 
chemical potential. The precise definition of the relevant 
’’single” electronic excitations £ a p (V) will be dependent 
on the specific model considered and is specified below. 
The index a of the spinors ip refers to valley/spin degrees 
of freedom. The partition function is 

Z = Tr e~ H/T = J Dip + Dipe~ s \' l ’«'‘ 4 ’ a \ ( 2 ) 

with measure defined by independent Grassmann vari¬ 
ables Dip + = n<A/j+, Dip = n dip a . The Matsubara ac- 

a a 

tion reads: 

r 1 / 1 ' r 

S [ip+,tp] = dr ip+ (r, r) (d T + K a p ) tpp (r, r) 

J 0 Jr 

-1 i’i (u r) ip+ (r, r) ipp (r, r) ip a (r, r), (3) 
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with the anti-periodic conditions, 

i>i(T,r) = -il)+(T + l/T,r), (4) 

ip a (r, r) = -ip a (r + 1/T, r). 

The local interaction term is not the most general one, 
but gen eralization to more ’’exotic” local cases (inter¬ 
valley 27 the exchange spin - spin couplin^H) i s quite 
straight forward. 

The normal and anomalous Green’s functions (consis¬ 
tent with definition in refPi) are: 

(r^ a (X) ip+ (X')) = ~G a0 (x-x ’); 

(X) V/3 (*')) = F a p (X-X '); (5) 

(t^+ (x)^+ (X')) = F+ 0 (X-,X'), 

where X = (r,r), X' = (r ',t'). For simplicity, we denote 
(i/j a (X)) as ip a (X), and drop the time ordering oper¬ 
ation T in correlators. For example, ( ip a (X) ipp (X')) 
stands for (Tip a (X) ipp (X')). Their expressions via par¬ 
tition function, Eq.([3]), are given in Appendix A. 

B. Space and spin rotations symmetries 

Generally a system may be invariant under both the 
space rotation and the spin rotation separately. Cer¬ 
tain valley symmetries are generally present. In Weyl 
semimetals, the action is typically only invariant under 
the combined space rotation and spin/valley rotation. 
The space rotation, r' = Ar, acts on a generalized spinor 
field as: 

V4 (r',r) = S (X) aa , x (r,T ), ( 6 ) 

^'p T ) = ( r ’ T ) S ( A )/3'/3 • 

The invariance of the action under the transformation, 
the correlators satisfy: 

G a p (X 1; X 2 ) = S (A) aa , G a ,f,, (A 1; A 2 ) 5 (A)J ; 

F aP (AG, X 2 ) = S (A) aa , F a ipt (A 1; A 2 ) 5 (A)^ (7) 

or in the matrix form: 

G (Ar, X 2 ) = S (A) G (A 1; A 2 ) 5 (A) f , ( 8 ) 

F (A lt A 2 ) = S (A) F (AG, A 2 ) 5 (A)* 

Assuming that the ground state is homogeneous, 
G(X, X) = G c and F (A, A) = p _1 A are constant ma¬ 
trices satisfying 

G c = S (A) G C S (A) + , A = 5 (A) AS (A)*. (9) 

Let E be generator of S (A). Then 

[E, G c ] = 0, EA + AS* = 0 (10) 


G c and A also satisfy the following equations: 

G+ = G C ,A‘ = -A. (11) 

In superconductor G c a0 = n c S a p and leads within the 
BCS approximation just to renormalization of the chem¬ 
ical potential. Therefore we finally obtain the Gor’kov 
equations, 

(■ ~d T S a p - K a p (V)) Gp 1 (A, A') -gx 
F 0a (A, A) F+ (A, A') = S (A - A') <5^; (12) 

(d T S a0 - K 0a (-■V)) F+ (A, A') - 
gF+ 0 ( A, A) Gp 1 (A, A') = 0. 

We will also discuss the non-superconducting state like 
the exciton condensate with opposite bulk properties, 
Gap (A, A) ^ n c 5 a p , F a p (A; A') = 0. In this case the 
Dyson - Schwinger form is more convenient. The gap 
equation can be recasted (see Appendix A) in matrix 
form as 

G-^G^+gSiX-X ')G' C , (13) 

where G' a0 is the traceless part of G a p, 

G' a p (A; A') = - (y a (A) Xp {X’)) c + (14) 

\s a pJ2(^ (X)^+(X')) c , 

with renormalized chemical potential taking case of the 
trace as in superconductor. 

The Matsubara Green’s functions (r is the Matsub- 
ara time) for uniform superconducting states can be ex¬ 
pressed via Fourier transforms, 

G a p (r, t; r',T r ) = T D [ d D kx 
(27rj J 

e ik-(r-r') ^ ^(r-F) . (15) 

n 

e ik.(r-r')^ e -^ n (r-r') F t^ (k)Wn) , 

n 

where w = ttT (2n + 1) is the Matsubara Fermionic fre¬ 
quency. The Matsubara Green’s functions in Fourier 
forms can simplify the calculation significantly. 

General formula 29 obtained in Appendix B for the en¬ 
ergy of the superconducting state reads: 

77r^ Tr{AA,) (16) 

while for chiral symmetry breaking (exciton condensa¬ 
tion) states the gaussian energy iP2> 

n = —Tr{— lnG + [Gq 1 G - 1]} + |l \G c G\ . (17) 
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III. THE DIRAC MODEL AND ITS 
SYMMETRIES 

A. Hamiltonian of the time reversal invariant TI 
with local interaction. 

Electrons on the surface of a TI perpendicular to z axis 
are described by a Pauli spinors ip (r), where the upper 
plane, r = {x,y}. In principle there are multiple valleys. 
The case of just one valley describing surface of the top o- 
logical insulator like Bi^Te^ were considered irf 1 13 SI It 
breaks time reversal invariance and does not allow chiral 
symmetry breaking, so here we consider the simplest case 
of multiple valleys: the Dirac model in which chiralities of 
the two Weyl modes are opposite described by field oper¬ 
ators ipf s (r), where f = L,R are the valley index (pseu¬ 
dospin) for the left/right chirality bands with spin projec¬ 
tions taking the values s =t,4- with respect to, for exam¬ 
ple, z axis. To use the Dirac (’’pseudo-relativistic”) nota¬ 
tions, these are combined into a four component bi-spinor 

creation operator, ip^ = ^ 14 , V 4 it> V’Jy.); whose in¬ 

dex 7 = {/, s} takes four values. The non-interacting 
massless Hamiltonian with Fermi velocity v F and linear 
dispersion relation, see Fig.l, readJ^S, 

£ 7 <5 = -ihv F V l al lS , (18) 

where two 4x4 matrices, i = x,y, 

“=(o -a)’ <“) 

are presented in the block form via Pauli matrices a. 
They are related to the Dirac 7 matrices (in the chiral 
representation, sometimes termed ’’spinor”) by a = /3 7 
with 



B. Continuous symmetries 

Symmetries of the 2D Dirac model with local interac¬ 
tions Eqs.fl l|18 ) were thoroughly discussed in relation to 
graphene 3 -^! They include parity, time reversal and the 
discrete chiral (flavor) transformation: ip — > 75 ip', ip —► 
—ip'Ys, where 75 = 70717273 = 75 - Spontaneous breaking 
of this symmetry has been comprehensively investigated 
in the context of graphen^H and will not be stressed here. 
There are also three continuous symmetries that in prin¬ 
ciple can lead to ordered phase with massless Goldstone 
bosons (order parameter waves). The first is the usual 
electric charge U (1), ip -» e lx ip, that is spontaneously 
broken in a superconducting state studied in next sec¬ 
tion. In addition there is the ’’chiral” flavour rotations 
SU ( 2 ) that play an important role in exciton condensa¬ 
tion that will be addressed in section V. 

1. Space time symmetries: just Aphelian space rotation 
combined with spin (pseudospin) rotation 

Unlike the non-interacting model the pseudorelativistic 
2+1 dimensional Lorentz invariance is explicitly broken 
by the static interaction Eq.([T]), so that only 2D rota¬ 
tions accompanied by the (pseudo) spin rotation already 
described in Section II with spin E operator, 

E =(o“,)=^' < 23 > 

remain a symmetry. The conserved quantity is the 
abelian angular momentum 

j = (r) jieij-uVj + j ip (r) . (24) 

The second part is referred to as the spin rotation, S = 

ISiP^ip. 

2. Electric charge 1/(1) 


The noninteracting Hamiltonian in these notations 
reads: 


The usual electric charge U (1) with conserved electric 
charge: 


K = ip t {-iHv F /3j ■ V — /i} ip. 


( 21 ) 


The Matsubara action, Eq. is conveniently written in 
pesudo-relativistic notations with ip = V’ + 7o with (Eu¬ 
clidean) 70 = i/3: 


Q = J ^ (r) ip (r) = ^ p (r). 


(25) 


Action Eq. (22) is invariant under the phase (global 
gauge) transformation, ip —> e lx ip 


s = [ f {-ip ('jod T + hv F 'Y- V-/x7o)VH- (22) 

J t —0 J r 

Since our focus is on symmetry and its spontaneous 
breaking, let us review known discrete and continuous 
symmetries. 


3. Flavour (chiral or valley) SU( 2) 

It was noticed early on in relation to graphen^23 that 
there is flavour SU (2) symmetry. It is shown in Ap¬ 
pendix C that quantities 

Qi = J* ^ ( r ) T i^ ( r ), (26) 
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FIG. 1. Schematic picture of the band reconstruction of Weyl semi-metal. 


commute with Hamiltonian and thus are conserved quan¬ 
tities.The generator matrices, 

Ti = -K,T 2 = - 75 , T 3 = - 7375 , (27) 

(75 = 70717273 ) constitute a nonrelativistic SU (2) alge¬ 
bra: 


[Ti,Tj] = i£ijkT k . 


(28) 


A discrete chiral symmetry is just the chiral rotation 
by angle 7r.The action Eq.(22) is invariant under in¬ 
finitesimal transformation, Sip = iTi%p\ Sip = iip'joTi'Yo- 
All three continuous symmetries commute?? (meaning 
charge, Flavour, and space time symmetry commute). 


C. Spontaneously broken electric charge U (1) 
symmetry phases: superconducting pairing channels 


Due to locality of the dominant interactions, the su¬ 
perconducting order parameter is local, 

0 = ( r ) Mapipl (r), (29) 

where the constant matrix M should be a 4 x 4 anti¬ 
symmetric matrix. Due to the rotation symmetry they 
transform covariantly under infinitesimal rotations gen¬ 
erated by the spin £ operator Eq. (231. 


Out of 16 matrices of the four dimensional Clifford al¬ 
gebra six are antisymmetric. We will not consider rather 
exotic phases in which in addition to the charge U (1) 
symmetry neither the 2D rotations or the SU (2) chi¬ 
ral transformations are spontaneously broken. Therefore 
superconducting order parameter is invariant under the 


remaining symmetries: [O, J] = [O, Qi ] = 0. Namely it 
is invariant under the SU (2) and is either scalar or pseu¬ 
doscalar under rotations. The requirement of invariance 
expressed using Eq. (23) takes a form: 


[O, J] = / ipt (r) M af) ipl (r) , 

J r,r' 

ip\ (r') £ ytjips (r')] (30) 

= - j ip 1 * (r) (EM + ME 4 ) ip ] (r) = 0. 

Similarly 

[O, Qi] = ~f ^ (r) {T t M + MT-) ^ (r) = 0 . (31) 


One finds that the only scalar is, M = ia y . There is also 
a pseudoscalar that will not be discussed here, namely 
we assume that the superconducting state preserves all 
the other symmetries. Which one of the condensates is 
realized at zero temperature is determined by the pa¬ 
rameters of the Hamiltonian along the line of dynamical 
calculation presented in section IV for the scalar. 

It turns out that there when the effective electron at¬ 
traction is replaced by repulsion and the superconductiv¬ 
ity is not realized there is still a possibility of continuous 
symmetry breaking that also belongs to a chiral univer¬ 
sality class: the chiral SU (2). 


D. Chiral SU (2) broken excitonic phases 

In this subsection we consider an opposite situation 
when the charge symmetry is unbroken, that is no su- 
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perconducting condensate appears. Still due to nontriv¬ 
ial multicomponent situation with the SU ( 2 ) symmetry 
there are possible transitions into a gapped exciton con¬ 
densate phases. It is plausible that rotational symmetry 
is also unbroken. Still there are two possible patterns, one 
is breaking down to an U (1) subgroup with two Gold- 
stone boson modes and another down to trivial subgroup 
with three Goldstone models. 

General order parameter now is 

p = J i’i ( r ) V a pil>f) 0 ), (32) 

where the constant matrix V should be an 4 x 4 hermitian 
matrix. There are four chiral SU ( 2 ) triplets of order 
parameters P = {P\,P 2l P 3 } (that can be viewed as the 
O (3) vectors). Their commutations with chiral rotations 
generators Qt defined in Eq.{26| are: 

[Pi, Qj\ = Jtf (r) [V h Tj] $ (r) = e ijk P k . (33) 

Four sets of matrices V a p in terms of the Dirac and chiral 
symmetry matrices are 

V (1) = {7o?2,7o7i, 2*^°} = 2 {-717273,*7o73,*7o} 

G (2) = {T 1 ,T 2 ,T 3 } = i {*73,7 o7i7273, -707172} 

P (3) = {71 ,7i73,7o7273}, (34) 

V {4) = {72,7273,707173} • 

There are also four chiral scalars, [ P , Qj] = 0, 
Q = (charge), S = 7172 ^ (spin), Hi = 

f = 7 o 7 i,* = x, y, that complete the Clifford 

algebra consisting of 16 hermitian matrices. These are 
not order parameters and hence will not be of interest 
to us. We also limit ourselves to the rotation invariant 
phases. The requirement of invariance expressed using 


Eq. (23) takes a form: 


[P,J} = t (r) . (35) 


Since V (3) and V (3) in Eq. (I34J) are not invariant under 
rotations only the first two are considered. If only one 
of the chiral vector order parameters has a nonzero ex¬ 
pectation value, say ^ 1 ' ^ 7 A 0 , the symmetry breaking 


pattern is SU (2) —► U (1), since 


p 3 (1) ,Q 3 


= 0. Accord¬ 


ing to the Goldstone theorem there are two soft modes in 
directions Q 1 and Q 2 . If in addition the second vector or¬ 
der parameter acquires VEV, ^P ? (2, ^> 7 A 0 and i 7 ^ 3, the 

pattern will be SU (2) —► / with three Goldstone modes. 

Which symmetry breaking mode is actually realized 
at given parameters of the system (chemical potential, 
Fermi velocity, interaction sign and strength, tempera¬ 
ture...) is a dynamical question. Therefore now we turn 
to dynamical aspects of the phase diagram of the Dirac 
model. 


IV. SUPERCONDUCTING STATE 

Within gaussian approximation, the Green’s functions 
obey the Gor’kov equations derived iii- 13 ' and in last sec¬ 
tion. For g > 0 anomalous Green’s functions are nonzero, 
so that the Gor’kov equations for Fourier components of 
the Greens functions simplify considerably, 

D~ 4 G Pk (w,p) - A 7/3 P* k ( U ,p) = <F K ; (36) 

D M F U + ^7/3 G /3« = °> 

where D~p = (ito — p) Syp—VFEijPiOt^p where the chem¬ 
ical potential is renormalized. The matrix gap function 
can be chosen as (A real) 

A^ = gF yP (0) = ( _° A 0 ) . (37) 

These equations are conveniently presented in matrix 
form (superscript t denotes transposed and I - the iden¬ 
tity matrix): 

D- 1 G-AF*=I; (38) 

P t - 1 P t + A*G = 0. 

Solving these equations one obtains 

G _1 = -D -1 + AD* A*; (39) 

P f = -_D 4 A*G, 

with the gap function found from the consistency condi¬ 
tion 

A* =-gJ2 Dt A*G. (40) 

ujq 

The off-diagonal component of this equation is: 

1 = ^ (a 2 + v 2 f p 2 + p 2 + h 2 uj 2 ) x 

9 top 

(a 2 + h 2 u} 2 + (vfP - X (41) 

^A 2 + H 2 oj 2 + (vfP + m) 2 ) 

The spectrum of elementary excitations obtained from 
the poles of the Greens function coincides with that found 
within the Bogoliubov - de Gennes approactf®: E p = 

±\JA 2 + (v F p- p) 2 . 


A. Zero temperature phase diagram for the 
superconductor - normal transition. 

At zero temperature the integrations over frequency 
and momentum limited by the UV cutoff A result in 

G=^A^-|log^|p4^, (42) 

z V A z + p z - p 
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FIG. 2. Order parameter at zero temperature as function 
of chemical potential of the TI surface Weyl semi-metal at 
various values of coupling parametrized by the renormalized 
energy U, Eq.(43l. For positive U (blue lines) the supercon¬ 
ductivity is strong and does not vanish even for zero chemical 
potential. There exists the critical coupling, U = 0 (red line), 
at which the second order transition occurs at quantum crit¬ 
ical point g = 0. For negative U the superconductivity still 
exists at g > 0, but is exponentially weak. 


where the dependence on the cutoff is incorporated in the 
renormalized coupling with dimension of energy defined 
as 


U = vf A 


47 rh 2 Vp 
9 


(43) 


This can be interpreted as an effective binding energy of 
the Cooper pair in the Weyl semi - metal. We consider 
only g > 0, since the particle - hole symmetry makes 
the opposite case of the hole doping, /i < 0, identical. Of 
course the superconducting solution exists only for g > 0. 
In Fig. 2 the dependence of the gap A as function of the 
chemical potential g is presented for different values of 
U. 

For an attractive coupling g stronger than the critical 
one, 


47 tH 2 vf 
9c = ’ 


(44) 


(when U > 0), blue lines in Fig. 2, there are two quali¬ 
tatively different cases. 

(i) . When g « U the dependence of A on the chem¬ 
ical potential is parabolic, se^ 13 l In particular, when 
g = 0, the gap equals U. As can be seen from Fig. 
2 , the chemical potential makes a very limited impact in 
the large portion of the phase diagram. 

(ii) For the attraction just stronger than critical, g > 
g c , namely for small positive U, the dependence becomes 
linear, see red line in Fig. 2, A = 0.663 g. So that the 
already weak condensate becomes sensitive to g. 

The case (i) is more interesting than (ii) since it ex¬ 
hibits stronger superconductivity (larger T c , see below). 
Finally for g < g c , namely negative U (green lines), the 


FIG. 3. Transition temperature as function of chemical po¬ 
tential at supercritical (U > 0, in blue), critical (U = 0, in 
red) and subcritical values of coupling. 


superconductivity is very weak with exponential depen¬ 
dence similar to the BCS one, A ss g exp[— (\U\ /g — 1)]. 
As was mentioned above, in the more interesting cases 
of large A the dependence on the chemical potential is 
very weak. A peculiarity of superconductivity in TI is 
that electrons (and holes) in Cooper pairs are created 
themselves by the pairing interaction rather than being 
present in the sample as free electrons. Therefore it is 
shown that it is possible to neglect the effect of weak 
doping and consider directly the g = 0 particle-hole sym¬ 
metric case. This point in parameter space is the QCP— 
and will be studied in detail in what follows. Of course, 
at finite temperature at any attraction, g > 0, there ex¬ 
ists a (classical) superconducting critical point at certain 
temperature T c that is calculated next. 


B. Dependence of the critical temperature T c on 
strength of pairing interaction. 

Summation over Matsubara frequency and integrations 
over momenta in the gap equation, Eq.( |42[ ), at finite tem¬ 
perature and arbitrary chemical potential. The critical 
temperature as a function of g and (positive) U is ob¬ 
tained numerically and presented in Fig. 3. Again at 
relatively large U the dependence of T c on the chemical 
potential is very weak and parabolic. When 0 < g < g c 
the critical temperature is exponentially small albeit 
nonzero, 0 « 1 + ({j) 2 - 

C. Zero chemical potential g = 0. 

At zero chemical potential the Hamiltonian Eq. 0 
possesses a particle - hole symmetry. Microscopically, 
Cooper pairs of both electrons and holes are formed, see 
Fig. la. The system is unique in this sense since the 
electron - hole symmetry is not spontaneously broken in 
both normal and superconducting phases. Supercurrent 
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FIG. 4. Phase diagram of STL Order parameter as function 
of U describing the deviation from criticality near the quan¬ 
tum critical point at A = 0, /r = 0. The critical line is a strait 
line in mean field approximation. 


in such a system does not carry momentum or mass. Per¬ 
forming the sum and integral over momenta in the gap 
equation, Eq.(42), analytically (see Appendix A), it be¬ 
comes (using the definition of U given in Eq.(43l) for 
U > 0: 


U = 2T log 


A ' 

2 cosh — 


(45) 


At zero temperature A = U, while A -> 0 as a power of 
the parameter U oc g — g c describing the deviation from 
quantum criticality 


T r = 


2 log 2 


U* 


zv = 1. 


(46) 


Here z is the dynamical critical exponent 18 . Therefore, 
as expected, the renormalized coupling describing the de¬ 
viation from the QCP is proportional to the temperature 
at which the created condensate disappears. 

The temperature dependence of the gap reads: 


A(T) = 2T cosh" 1 Q exp . (47) 

This it typical for chiral universality classed®. 

It is interesting to compare this dependence with the 
conventional BCS 291 for transition at finite temperature, 
namely away from QCP, see Fig lb. At zero temper¬ 
ature A (0 )/T c = 21og2 = 1.39 (within BCS - 1.76), 
while near T c one gets A /T c = 2 3 / 2 log 1 / 2 2\/l — t = 2. 


35-/T ~t (BCS - 3.07vT^I), where t = T/T c . To de¬ 
scribe the behavior of the STI in inhomogeneous situa¬ 
tions like the external magnetic field, boundaries, impuri¬ 
ties or junction with metals or other superconductors, it 
is necessary to derive the effective theory in terms of the 
order parameter A (r), where r varies on the mesoscopic 
scale. 

Phase diagram of STI order parameter as function of 
U at zero temperature is plotted in Fig.4. 


D. Coherence length and the condensation energy 


The quadratic term of the Ginzburg-Landau energy 
F-2 = A*r (p) A p is obtained exactly from expanding 
the gap equation to linear terms in A for arbitrary exter¬ 
nal momentum. The dependence on p is non-analytic 
and within our approximation higher powers of p do 
not appear. The second term is very different from the 
quadratic term in the GL functional for conventional 
phase transitions at finite temperature 19 or even quan¬ 
tum phase transitions in models without Weyl fermions 18 
and has a number of qualitative consequences. Compar¬ 
ing the two terms in Eq.(51), one obtains the coherence 


length as a power of parameter U oc g — g c describing the 
deviation from criticality: 


€(U) = -v F KU- 


v=\. 


(48) 


This is different from the dependence in non-clriral uni¬ 
versality classes that iP^£ (T) oo (T c — T)~ v , v = 1/2 in 
mean field. Of course in the regime of critical fluctua¬ 
tions this exponent is corrected in both non-chiraP 9 and 
chiraPS universality classes. 

Local terms in the GL energy density are also calcula¬ 
ble exactly. Expression for the kernel can be written as 
a trace: 



1 ^ h 2uj2 — v fP ' 9 + 

9 u, q (h 2 uj 2 + v 2 F q 2 ) (h 2 uj 2 + v 2 f |q-p| 2 ) 


(49) 


Integrating over u> (at zero temperature) one obtains 


1 I" pq cos (f> — p 2 
8ir 2 h 2 v F J q ^ | q - p| 2 + q |q - p| 
U 

Airh 2 v F ’ 


(50) 


where cj> is an angle between q and p. The integral is 
homogeneous in momentum and therefore is linear in p = 
|p| and one arrives at: 


r(p) 


u , \p\ 

4,ttH 2 v f 16 vfK 2 


(51) 
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E. Local terms in Ginzburg - Landau equation 
and energy 


For p = 0 the gap equation Eq. (42) reads 


Vk*(~ U+V A^)= 0 , 

This is obtained from the energy functional 

F =«kf/ <i2r {-' 74 ' a+ !( a-a ) 3/2 }- 


(52) 


(53) 


It is quite nonstandard compared to customary quartic 
term (A*A ) 2 in conventional universality classes. The 
GL equations in the homogeneous case for the conden¬ 
sate gives A 0 = UP with critical exponent 0 = 1, dif¬ 
ferent from the mean field value 0 = 1/2 for the U (1) 
universality class 19 . The condensation energy density 
is fo = — 127r n 2 v' 2 & 2 ~ a with a = —1. The free en¬ 
ergy critical exponent at QCP therefore is also differ¬ 
ent from the classical a = 0. The magnetic field cou¬ 
ples to the order parameter field by a standard minimal 
substitution. The Ginzburg-Landau approach the only 
available practical tool to study properties of inhomoge¬ 
neous configuration ^ 3 in external magnetic field like the 
Abrikosov vortex systems. 


V. EXCITON CONDENSATION 

A. Chiral symmetry breaking 

Assuming that the electric charge U (1) and rotation 
symmetry is unbroken (no superconductivity) for insuffi¬ 
ciently strong attractive interaction (pairing) one still can 
have transitions due to always existing effective repulsion. 
This possibility was already considered in graphene 32 . 
Here we note that the symmetry patten might be very dif¬ 
ferent from often invoked relativistic 2 + 1 Gross - Neveu 
model thoroughly studied as a toy model in relativistic 
quantum field theory. Possible chiral symmetry breaking 
states were reviewed in subsection lid and we consider 
first a ground state with nonzero order parameter 

(P 3 (1 ) )=u^0. (54) 

It is convenient to parametrize the gaussian variational 
ground state by the trace of propagator: 

G(X-X) = m/g 7o , (55) 

with ’’mass” m determining the order parameter. Ac¬ 
cording to gap equation derived in Appendix B, we have 

= ( ioj + p) I — a ■ ( HvFk ) — m+o; (56) 
£ _ (p + ioj) I + a ■ (, hvFk ) + myo 
m 2 + (hvFk ) 2 — (p + i oj ) 2 


leads to 


ml = gG' (t, r; r, r) 7 0 1 = - x (57) 


_ 9 _ 

( 2 ^r 


m 


Jui ,k m 2 + (hvFk) 2 - (p + i w ) 2 

Performing integrations with momentum cutoff A one ob¬ 
tains for m/ 0 : 

, _ 9 f (V m 2 + A 2 - m) , p < m , . 

47 t (hvp) 2 1 (Vm 2 + A 2 - p) ,p>m 

For p = to one obtains the first order transition point 
(for p > 0 ) 


4-7T ( Hvf ) 2 
\9c\ 


= \fpV+ A 2 - , 


(59) 


For larger repulsion, \g\ ^ \g c \, the mass is larger: 


1 ... ... /A 2 2 tt (hv F ) 4 \ 

m = 9 +~—-2 (W - Iffcl) I + i„i i_ / 1 • ( 60 ) 


(hv F ) 


8tt 


I g\ I Sc 


The free energy, using the Abrikosov formula derived 
in Appendix B, is (A » m) 


r m dV 

(fl (m) — fl 0 (m = 0)) /V = / 2m 2 ——dm (61) 

J o 


m 


I o 2tt(Hvf) 


r dm = — - 


67t ( Hvf)~ 


For p > to, since 


d(i/\g\) _ 

dm 


> 0 , the 


Ait(Jivf ) 2 \/m 2 + A 2 

chiral symmetry breaking state will have higher energy 
than the normal one. For \g\ < \g c \, the stable ground 
state is the normal state with m = 0 . 

Finite temperature properties including the phase di¬ 
agram and the QCP at zero chemical potential can be 
studied along the lines similar to the superconducting 
transition. 


VI. DISCUSSION AND CONCLUSIONS 

Having considered both the superconducting and the 
excitonic transitions for sufficiently strong attraction or 
repulsion within the gaussian approximation, a natural 
question is what happens when we approach the quan¬ 
tum critical point. At criticality various renormalization 
group methods should be usecpA 


A. Criticality beyond the gaussian approximation 

Critical (quantum) fluctuations are expected to be sig¬ 
nificant in this relatively low dimensional (relativistic 
2+1 dimensional system) system. Generally they are not 
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as strong as in 2D statistical system at finite tempera¬ 
ture, but stronger than in 3D one. The approximation 
we have made describes reasonably well ’’gaussian” fluc¬ 
tuation beyond the region where stronger critical fluc¬ 
tuations in these systems appear and should be treated 
nonperturbatively 20 typically using variants of the renor¬ 
malization group approach 19 . The critical exponents in 
this region differ from the one called ” quantum gaussian 
(BCS)” in ref.® and available results are obtained using 
either e expansioi l 21 * 34 ! (e = 4 — d, where d = 2 + 1 is the 
space-time dimension), 1/N, where N is the number of 
fermionic species on the surfac ^ 21 and functional (strong 
coupling) RG 35 and Monte Carlo simulation ^ 37 (with 
reservations specified below). The universality class of 
the supconducting transition according to classification 
proposed in ref.^U is the chiral XY (symmetry of order 
parameter U (1)) with N = 1. The large N expansion 
is not reliable for the one component system considered 
here (but the number might be larger in similar systems 
for which our approach trivially generalizes), so let us use 
the £ expansion. 

Using the formulas for the anomalous dimensions of 
the order parameter (see second reference irPU), 

1A = V = l/4er + 0.044s 2 + O (s 3 ) « 0.294 (62) 

and its square, 

7a2 = (l + \/li) /10s+0.065e 2 ~ 0.43s+0.065e 2 , (63) 

critical exponents are obtained from the hyperscaling re¬ 
lations: 

a = 2 - d/ (2 - 7 A 2 ) = (e - 2 7a2 ) / (2 - 7a2 ) = -0.353, 
P= (1 + 7a)/ 2(2- 7a2 ) =0.515. (64) 

and can be compared with those in Table 1 in ref® The 
exponents from the £ expansion were found to be consis¬ 
tent for larger values of N with the latest Monte Carlo 
simulation^!, while consistent with the functional RGP. 
The critical exponents of the chiral transition belongs to 
this class with N = 2. 

The corresponding chiral universality class for the ex- 
citon condensation in Dirac semi-metal is the Heisenberg 
N = 2 (SU (2)). The c ritical e xponents for this case 
were also calculated in ref J 21 l 34 l d S] Recently they were in¬ 
voked in a discussion of second order quantum transitions 
in Hubard model on honeycomb lattice®. It should be 
noted that the flavour symmetries are often broken ” ex¬ 
plicitly” by some kind of anisotropy. In this case Gold- 
stone bosons acquire a small mass (like pions in quantum 
chromodynamics in which the chiral symmetry is slightly 
broken by the light quark masses), although their major 
properties remain intact. This can be taken into account 
as a small perturbatiorP®. 


B. Experimental feasibility of observation of 
quantum phase transition 

The best candidate to observe the superconductivity 
is a topological insulator of the Bi 2 Se 3 family. To esti¬ 
mate the pairing efficiency due to phonons, one should 
rely on recent studies of surface phonons in Tl 12 . The 
coupling constant in the Hamiltonian, Eq.(Jl]), is obtained 
from the exchange of acoustic (Rayleigh) surface phonons 
g = Xvph 2 /27r/i, where A is the dimensionless effec¬ 
tive electron - electron interaction constant of order 0.1 
(somewhat lower values are obtained in ref.®). It was 
shown in ref.® that at zero temperature the ratio of A 
and /-I is constant with well defined g, —> 0 limit with value 
g = 0.23 eV nm 2 for vf ~ 7 • 10 5 ?n/s (for Bi 2 Se 3 ). The 
critical coupling constant g c , Eq.(|44|), can be estimated 
from the Debye cutoff T D = 200 K determining the mo¬ 
mentum cutoff A = Td/c s , where c s is the sound velocity. 
Taking value to be c s = 2 • 10 3 m/s (for Bi 2 Se -$), one ob¬ 
tains g c = Attvfc s H 2 /T d = 0.20 eV nm 2 . 

Of course the Coulomb repulsion might weaken or even 
overpower the effect of the attraction due to phonons, so 
that superconductivity does not occur. In TI like Bi 2 Ses 
however, the dielectric constant is very large e = 50, so 
that the Coulomb repulsion is weak. Moreover it was 
found in graphene (that has identical Coulomb interac¬ 
tion) , that although the semi-metal does not screen®, the 
effects of the Coulomb coupling are surprisingly small, 
even in leading order in perturbation theory. The super¬ 
conductivity was observed in these systems that howeved 
had to be either doped in the bulk or on the surface® 
(by a Cu) or by applying pressure®. It is not yet clear 
whether the observed superconductivity is a bulk or a 
surface effect. 

The Dirac semimetal in optically trapped cold atomd® 
offers a well controllable system in which this phe¬ 
nomenon occurs both for repulsive interaction (chiral 
symmetry breaking) and in particular the attractive one 
(superconductivity) becuase there is no Coulomb repul¬ 
sion as the atoms are neutral. 

Recently after experimental discovery of 3D Dirac 
semi-metald® the new class of questions similar to 
those discussed in present paper arise. Extraordinary 
electronic properties of these Dirac material^® includ¬ 
ing superconductivity® and chiral condensate are being 
studied theoretically and experimentally. 


C. Conclusions 

We have studied continuous phase transitions in a 
Dirac semi-metal realized recently as a surface of topo¬ 
logical insulator. The noninteracting system is character¬ 
ized by (nearly) zero density of states on the 2D Fermi 
manifold. It degenerates into a point when the chemi¬ 
cal potential coincides with the Weyl point of the sur¬ 
face states as in the original proposal for a major class 
of such material^®. The pairing attraction (the most 
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plausible candidate being surface phonons) therefore has 
two tasks in order to create the superconducting conden¬ 
sate. The first is to create a pair of electrons (that in the 
present circumstances means creating two holes as well) 
and the second is to pair them. To create the charges 
does not cost much energy since the spectrum of the Weyl 
semimetal is gapless (massless relativistic fermions); this 
is effective as long as the coupling g is larger than the crit¬ 
ical g Cl see Eq.(44). The situation is more reminiscent of 
the creation of the chiral condensate in relativistic mass¬ 
less four - fermion theory (a 2 D version 2 ^ was recently 
contemplated for graphen tffl -— than to the BCS or even 
BEC in condensed matter systems with parabolic disper¬ 
sion law. Due to the special ’’ultra-relativistic” nature of 
the pairing transition at zero temperature as a function 
of parameters like the pairing interaction strength is un¬ 
usual: even the mean field critical exponents are different 
from the standard ones that generally belong to the U (1) 
class of second order phase transitions. 


To summarize, we studied the phase diagram of the su¬ 
perconducting and chiral transition at arbitrary chemical 
potential, effective local interaction strength and temper¬ 
ature T. The quantum (T = 0) critical point appears 
at zero chemical potential and belongs the Un (1) chiral 
universality class (the subscript denotes number of mass¬ 
less fermions at QCP according to classification ii i — I— () 
for the attraction (superconductivity) and SUn (2) for 
repulsion (exciton condensation). 
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Appendix A: Path integral derivation of the Gorkov 
and Dyson - Schwinger equations 


1. General correlations and sources 


We introduce the grassmanian source terms into Mat- 
subara action Eq.([3|: 


[</’+,v>] 


(Al) 


J dr j [4>+ (r, r) J a (r, r) + i/j a (r, r) J+ (r, r)] . 


The generating functional for the disconnected correla¬ 
tions is Z (J + , J) 


(V’a (X)) 
(A)^ (X')...) U 


1 


6Z 


Z [-5J+ (A)] ’ 

1 5Z 

% [-SJ+ (A)] \-SJ+ (A' 


(A2) 


where A = (r,r), A' = (r ' ,t'). The correlators is (Mat- 
subara) time ordered. The generating functional for the 
connected correlations are W = — In Z (J + , J): 


PO) 

<Ttfi a (A)^ (X')...) 


SW 

sj+Jx)’ 


(A3) 


[-SJ+ (A)] 


SW 

~-5J+ (A') 


In particular, we define the normal and anomalous 
Green’s functions (see refP*) as, 

/ , \ SW 

(t*‘ (x > a m ) = {XX); 

SW 

x »= - sjtwsjt rn = {XX) • 

(tv: m # <*'>)=- < X;X ' iA4 > 

where we denote (ip a (X)) as ip a (X), and drop the 
time ordering operation T in correlators. For example, 

(ipa (A) ipp (AO) stands for (Ttp a (A) ipp (AO). 


2. Derivation of the Gor’kov equation for 
superconductivity 


Using identities, 

< A5 > 

the equations of motions lead to 

(hd T S a p + K aP (V)) 1/ip (A) g x 
(typ (A) fPp (A) W (A)) = J a (r, r) , (A6) 

(tidrSaf, - kp a (-V)) ^ (A) +gx 
(</£ (X)^(X)^ (A)) = J+(r,r). 

Generally full correlations of the fields can be expressed 
via connected correlators. For example, 
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(X) 4> p (X) </> Q (X)} = ^ (X) ^ (X) v>« (X) - 

1>a PO PO P0) c + ^ PO X (A7) 

(X) (X)) c + (X)) (X) ^ (X)) 

+ (X) ipp (X) ip a (X)) c . 

For Gaussian mean field approximation, we 
omit higher order connected correlations, like 

(X)^p (X)^ a (X)), 

J a (X) = (ficU a/ 3 + K a p) 1>p (X) -g x (A 8 ) 

1 # (X) ^ (X) Vg PO - ^ (X) (X) P0) c 

+# (X) (</, a (X) ^ (X)) c +(X) (X) Vg (X)) c 

Performing functional derivative SJ j x ^ (by using the 
identity, gjjx^ ^a (X) = (ip a (X) ip+ (X')), the corre¬ 
lators of odd number of Grassmannians vanishing) and 
taking at the end J a (X) = 0, one obtains: 

s ( X-X ') <s Q7 = (s T< y a/S + x a/3 ) (X) (X')> -9 

x {- (X) V>+ (XO) (X) (X)) (A9) 

+ <^ (X)^+ (x')>(vg P0tf£ PO) 

+ (X) V-+ (x')) (X) VG (X))}, 

or 

<5 (X - X') <5 a7 = (~d T S a p - K a p} Gp 1 (X, X') -g 

x {—G a7 (X, X') Gpp (X, X) (A10) 

+G /37 (X, X 7 ) G aj g (X, X) c + 

F+ (X,X')F> a (X,X)}. 

Similarly the second equation of motion, 

Ja ( T > r ) = (d T Sat 3 - Kp a (-V)) (t, r)^ +g 

{{^a (. T x)ipp (T,r)^> (ipp(T,r)) (All) 

- ( T > r )) (V’/S (t, r ) (T, r)^ 

+ (T,r))(^(T,r)^ ( T > r ))|, 

gives 

(-V)) F+ (X, X') + 5 {-E+ (X, X) 
G/j 7 (X, X') + F + 7 (X, X') Gpp (X, X) + (A12) 

Gp a (X,X)F+ (X,X')} = 0 . 


In the superconducting phase, if the spin rotation in¬ 
variance is not broken and chiral symmetry is preserved, 
G a p (X,X) = G c a p = n c S a p. For the N component 
spinors n c = n/N, where n is the density of electrons. 
The quadratic parts of Eqs.( A10|A12 ) simplify: 


-g {-G ai (X, X’) Gpp (X, X) + Gp 1 (X, X') G a p (X, X)} 
= -9 {~G ai (X, X') Nn c + n c G ai (X, X')} 

= 5 (X-l)n c G« 7 (X,X , ) ; (A13) 

g {X+ (X, X') Gpp (X, X) + Gp a (X, X) F+ (X, X')} 

= -g{N-l)n c F+(X,X'), 


and such terms can be absorbed to the chemical potential 
term with the chemical potential replaced by the renor¬ 
malized 


g + g(N- l)n c = Hr. (A14) 

Therefore we finally obtain the Gor’kov equations, given 
in Eqs.(fl2l). 


3. Derivation of DS equation and renormalized 
chemical potential for nonsuperconducting state. 

We need only to discuss the equation for J^j-, as |^T is 
just the complex conjugate of and ^ ^ = 0 . 

We reorganize the equation f^J-as 

= ( d + K a p) 5 (X - X') - g5 (X - X') x 

oipp 

| l^ a (X) (X)) c - \& a p ^ (X) (X)) c j 

+g6 a p6 (X - X') {fa (X) i,+ (X)) c - (A15) 

£<^P0^P0) c }> 

The last two terms are proportional to 5 a p , and can 
be absorbed to the chemical potential, 

H T 9&a/3 {£ (V> 7 (X)^(X)) C (A16) 

- (if>p (X)^+ (x)) c }= M «. 

The gap equation can be recasted as 

G -1 = Gq 1 + gS (X - X') G' (X; X') 

Go 1 = ( d T S a p + K a p ) <5 (X - X') (A17) 

Gap (X; X') = - { (^ a (X) $+ (X')) c - 

^^(^(Xl^KX'))), 

where G( 1;/ , ; is the traceless part of G a p , and the chemical 
potential in K a p is hr. 
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Appendix B: A formula for energy of the 
superconducting state 

1. Gap equation in the Nambu notation 


A compact representation of the Gorkov equations for 
superconductors is the Nambu notations 


Ga0 (X,x') 


so that 


/ Sj> a (X) S^ a (X) \ 
1 5 Jfi(X’) S J+( A'') \ 
W+(X) S^j(X) 
\SJp(X') 6J+(X')J 

(~G a 0 

v F ip(. x i T i x 'i T ') 


(Bl) 

F a 0 {x,T\ x' ,t')\ 
Gf 3 a (x',t';x,t)J 


g-i-f&o 1 ) u 0 

(Go 1 )!! = (d T 5 a p + K a0 (V)) 5 (x — x' ,t — t’) ; (B7) 
(Gq 1 ) 22 = (d T Sa0 ~ Kp a (-V)) S (x - x' ,T - t') , 

and the gap equation can be cast in the Dyson form 


0 S 12 

S 21 0 


= Go 

S 12 = gS (x — x', t - t') (ip a (r, r) (r, r)) c ; (B8) 
S 2 i =gS(x - x', t - t ') (t/)+ ( t , r) (r, r)^>. 


( XLsl 

I Wp 

I SJ i 

\Wp 



(B2) 


2. Derivation of the expression for energy density 

The free energy is 


.The functional identity 


SJ a (X) <% SJ a (X) 

di/jp SJ 1 ( X') 8 ipp 


SJ-y (X') 


8 ai (X - X '), 


where abbreviations 


(B3) 


SJ a (X) 5^0 _ v r SJ a (x) 

8J 7 (X>)-^J ' 5^0 (X") 


(B4) 


(X") 
SJ. 7 (X') ’ 


Say (X 


X') = 5 ai 5 {X - X ’), 


are used. In the Nambu matrix form it reads 


8 J„ SJ. \ / 8i/>p Stjjp \ 

We Wf \ I Jx; JJ+ \ 

SJ+ 5J+ I I 81 !>+ 5i/j+ 

We 81 P+J \T77 JJ+/ 


(Say 0 \ 

V 0 Say) ' 


(B5) 


The derivatives are: 


m = ~- p { ^Tr{- ln<? + [Q^Q - 1]} + <f>[G] 
where in Gaussian approximation, 



m = \ J F a 0 {X-X)F+ a {X-X). (BIO) 

The gap equation can be also obtained by X il[Q] =0. 

The energy difference between the superconducting 
state and normal state can be obtained by the differ¬ 
entiating of the grand canonical potential with respect 
to the coupling constant: The green function Q is 

dependent on g, but due to sg ^ [G, 9] — 0 , 

±n[G,g} = -^n[G,g] (Bll) 

= -^ j ^ F a 0 {X; X) F^ a (X; X). 

For homogeneous state, F a 0 (A; X ), Ft (A; A) are con¬ 
stant (not dependent on x,t), and we introduce 


+ Ka0) s (A - A'); 

5 J + =9S{ A - AO (^ Q (r, r) ip p (r, r)) ; 

(X') 

( x'] =gS(x - X ' } (T ’ r) ^ (T ’ r) ) 5 (B6) 

= & 6 - K W (-V)) <5 (A - AO • 

For the non-interacting model 


F a 0 (A; A) = ~A a 0,F+ 0 (A; A) 

= 1a+ = 1a* 

g a? g A W’ 


we obtain the free energy formula 33 

|n| 5 ]=- ^tv(aa+) 


dn[G]=^d Q) Tr(AA+) 


(B12) 


(B 13 ) 
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Appendix C: Gap equation and free energy for 
chiral symmetry breaking states 


2. Details of the calculation of the condensate and 
energy 


1. Derivation of the gap equation 

We will discuss the non-superconducting state with 
G a p (X,X) n c 5 a /s, F a p (X; X') = 0, which happens 
for example in the case of chiral symmetry breaking 
(CSB) state. We need only to discuss the equation for 

as^y^ is just the complex conjugate of and 


Using the integral 


SJ+ _ SJ a 

SJ, 


= 0. We reorganize the equation fjy^as 


5tpp 


= (d T 5 a p + K a p) S (X - X') — gS (X — X’) x 


/ ^ a (X)^+(X)) c -^5 a pX (Cl) 

E <^7 ( A ) E ( A ')) c } + 9 SapS (X - X 1 ) x 
{(# (X) ^ (X)) c - E (V’t P0 < P0> c }, 

The last two terms are proportional to 6 a p, and can 
be absorbed to the chemical potential, 

“Ws: (^(X)i,+ (X)) c - (C2) 

(V>/3 ( A '))J = PR- 

The gap equation can be recasted as 


G _1 = Gq 1 + gS (X — X') G’ (X; X') 

7-1 


G 0 1 — (yd T 5 a p + K a p^ SS (X — X') 

G' a p (X; X') = - { (y a (X) (X')} c - (C3) 

^^E(^( A )^( A ')> C 


where G' a p is the traceless part of G a p, and the chemical 
potential in K a p is /ifl.The free energy is now 

m = -(/?)" 1 {Tr{-lnG+[G 0 - 1 G-l]} (C4) 

+*m, 

where 

$[£] = -§ J G ' a p (x,T-,x,T)G'p a (x,T;x,T). (C5) 

A similar free energy equation can be obtained, 


|n[ S , 9 ] = |n[ S . 9 ] 


(C6) 


1 

2 P 


/ G 'ap ( x i T i x i T ) G 'pa i x i T i x i T ) ■ 
J r,r 


m 2 + (hvFk ) 2 — (/.i + * w ) 2 


(C7) 


:© ^ m2 + {hVpk ) 2 — 


and 


(flVFk)' 


(27r)‘ to 2 + 

2 yA p© (Vm 2 +p 2 - ^ 

4tt J p= q y' to 2 + p 2 


(C8) 


J (\/rra 2 + A 2 — rra) , // < m 

I 4 ^ (pTO 2 + A 2 - p),p>m 

the gap equation becomes: 

m = ff / (v to 2 +A 2 -to),/x<to . . 

47r(fiu F ) 2 1 (Vm 2 + A 2 - fj) , /j, > m 

Appendix D: Symmetries of the Dirac model 
1. The chiral nonrelativistic SU (2). 

Generally the charge algebra is [Q, Qi] = 0; [Qi, Q :) ] = 
iEijkQk- Indeed all three chiral charges commute with 
density: 

\P (x ), Qi] = [ [E (x) i/> (x), ^ (y) Till) (y)] (Dl) 
= ^ ( X ) [I,Ti]ip(x) = 0, 
and the kinetic term, 

[ H , Qi] = [ (x) (aipi + a 2 p 2 )i> (a), ^ (y) (y)] 


= V’ t [x] [aipi + a 2 p 2 , Ti] if) (x) = 0. 


(D2) 


The three generator matrices, Eq. (27), constituting the 


SU (2), [Ti,Tj] = iSijkTk, commute with both 7071 and 
7072 that appear in noninteracting Hamiltonian, Eq. ([Tj) . 
The density - density interactions part of Hamiltonian 
also commute with the SU (2) charges since 


[p(r),Qi]= [E {r)ip(r) ,^ f (r 1 ) T^ (r')\ (D3) 

J r' 

= ^ (r) [J, Ti] if) (r) = 0. 
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Correspondingly the action Eq.(22) is invariant under 


Sip = iTiip-, dip = iipjoTi'yo. Indeed both the energy 
term, 


Mn = - 


*V’7oT i 7o7oi9o'0 + ip S r yodoiT i ip (D4) 


-iipryoTidoi’ 8 + ip 'yod 0 iT i i/j s 


= 0 , 


and the momentum terms, 


5 Ai = - j [iip'yoTi'yo (7 • d) ip s + ip (7 • d ) iTiip] 

= -i J ip [70^707 + 7 TP\ ■ dip = 0, (D 5 ) 

are invariant. The nonrelativistic interaction term 

SA = v J [iipjoTi'yo'yoip + ipioiTiip] (ip r 7oV’ r ) (D6) 
= iv J [—ip'yoTiip + ip'yoTi'ip] (Vlo ip r ^j = 0. 
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